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ABSTRACT 

The  concept  of  item  by  item  sequential  sampling  has  been  developed 
from  the  Wald  sequential  probability  ratio  (SPR)  test.  The  SPR  plan 
for  one-sided  variables  inspection' (a  plan  fcr  testing  the  mean  X'  of  a 

— - — O 

normal  distribution  with  known  variance  o'  )  is  characterized  by  two 
parallel  lines  on  a  plot  of  the  cumulative  sum  of  observations  T(n) 
versus  accumulated  sample  size  n.  The  operation  of  the  plan  is  described 
by  analyzing  the  sequence  T(n)  For  the  case  where  small  values  of 
X'  are  preferred,  T(l)  =  X^  is  the  first  observation.  If  T(l)  falls  on 
or  above  an  upper  bound  Ta(l),  the  lot  is  rejected,  if  T(l)  falls  on  or 
below  a  lower  bound  T^l),  the  lot  is  accepted.  If  T^(l)  <  T(l)  <  Tu(  1 ) , 
another  sample  is  taken.  T(2)  =  X-j  +  X2  is  then  treated  in  a  similar 
manner,  and  so  on  until  a  course  of  action  is  defined. 

The  number  of  observations,  n,  required  for  defining  a  course  of 
action  is  a  random  variable,  often  characterized  by  large  values  of  n 
before  a  decision  is  made.  Truncation  rules  are  required  to  prevent 
unacceptable  sample  sizes  even  though  convergence  is  guaranteed.  This 
truncation  increases  the  associated  producer's  and  consumer's  risk. 

Herein  a  sequential  non-probability  ratio  (SNPR)  plan  is  examined 
which  is  characterized  by  intersecting  lines  and  thus  an  implicit 
truncation  property.  A  heuristic  algorithm  to  define  the  parameters  of 
the  plan  is  introduced.  The  algorithm  produces  an  SNPR  plan  which 
provides  the  same  levels  of  protection  as  the  non-truncated  SPR  plan. 
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CHAPTER  I 


INTRODUCTION 


Problem  Statement 

Abraham  Wald  developed  a  sequential  probability  ratio  (SPR)  test 
in  1943  to  test  statistical  hypotheses.  Initially  used  in  sequential 
acceptance  sampling  by  attributes,  it  has  since  been  adapted  to 
variables  inspection.  A  unique  feature  of  Wald's  SPR  test  is  that 
the  number  of  observations  required  by  the  procedure  is  not  determined 
in  advance  of  the  experiment.  Its  ability  to  quickly  detect  good  or 
bad  quality  lots  makes  it  superior  to  other  acceptance  sampling  pro¬ 
cedures.  It  can  materially  reduce  the  required  amount  of  inspection 
by  50%  on  the  average  when  compared  to  the  single  sampling  plan. 

While  it  has  been  shown  that  the  probability  equals  1  that  the 
procedure  eventually  terminates  [9],  there  is  no  effective  upper  limit 
on  the  number  of  items  to  be  inspected.  If  lots  are  mediocre  in  quality, 
the  number  of  units  required  to  be  inspected  can  be  quite  large.  This 
problem  is  compensated  for  through  the  use  of  truncation  rules  which 
prohibit  sampling  beyond  a  certain  number  of  units  by  forcing  a  decision 
at  some  a-priori  fixed  integer.  If  the  procedure  is  truncated  prema¬ 
turely,  its  ability  to  discriminate  between  good  and  bad  lots  is  weakened. 
If  the  procedure  is  truncated  too  late,  its  ability  to  discriminate  is 
only  marginally  reduced  but  it  loses  its  advantage  of  reducing  the  number 
required  for  inspection. 

This  study  provides,  through  a  procedure  proposed  by  Guild  [5], 
a  sequential  non-probability  ratio  (SNPR)  test  which  is  bounded  for 
the  case  of  variables  inspection.  The  specific  objectives  of  the 


thesis  are: 


1.  Development  of  a  heuristic  procedure  which  specifies  the 
arguments  of  an  SNPR  procedure  fQr  a  given  set  of  parameters. 

2.  Comparison  of  SPR  and  SNPR  procedures. 

3.  Development  of  a  computer  code  for  the  heuristic  and  SNPR 
procedure. 

4.  Generation  of  a  select  set  of  SNPR  test  plans. 


Introductory  Remarks 


The  purpose  of  an  acceptance  sampling  plan  is  to  define  a  course 
of  action;  i.e.,  acceptance  or  rejection  of  a  particular  lot  of  given 
quality.  If  all  lots  are  of  the  same  quality,  the  plan  will  indicate 
acceptance  of  some  lots  and  rejection  of  others,  and  the  accepted  lots 
will  be  no  better  than  the  rejected  ones.  If  the  lots  differ  in 
quality,  the  plan  will  accept  good  lots  more  frequently  than  it  will 
bad  lots  II]. 

Usually  a  plan  is  designed  so  that  material  considered  to  be  of 
good  quality  will  have  a  low  probability  of  being  rejected  (producer's 
risk)  while  material  considered  to  be  of  bad  quality  will  have  a  low 
probability  of  being  accepted  (consumer's  risk).  The  grade  of  material 
considered  good  is  called  the  "Acceptable  Quality  Level"  (AQL)  and  the 
probability  of  rejecting  material  of  grade  AQL  or  better  is  designated 
by  a.  Similarly,  the  grade  of  material  below  which  rejection  takes 
place  is  referred  to  as  the  "Lot  Tolerence  Percent  Defective"  (LTPD) 
and  the  probability  of  accepting  material  of  grade  LTPD  is  designated  as  8 


3 


The  precise  relationship  between  lot  quality  and  the  probability 
of  acceptance  (P  )  of  any  given  lot  is  shown  by  the  operating  charac¬ 
teristic  (OC)  curve  of  the  particular  sampling  plan.  The  curve  graphi¬ 
cally  illustrates  the  plan's  ability  to  discriminate  between  different 


quality  lots  with  fixed  o'.  The  pairs  (1-a,  X'^)  and  (b,  X'^p^), 
where  X ' 1S  the  mean  locating  a  percent  defective  equal  to  a  desired 
AQL  and  X'LTpD  is  the  mean  locating  a  percent  defective  to  a  desired 
lot  tolerance  LTPD,  specify  two  points  on  the  curve  as  shown  in 


Figure  1.  The  point  on  the  curve  where  P,  equals  .5  locates  a  mean 

a 

known  as  the  indifference  point  ( X '  ) .  At  this  point,  there  is  no 

a 

particular  preference  for  acceptance  or  rejection  of  the  lot  (a=e). 


The  curve  shown  in  Figure  1  describes  a  plan  where  low  values  of  X' 
are  desirable.  Besides  showing  lot  quality  relationships,  the  OC  curve 
is  also  useful  in  evaluating  alternate  sampling  plans. 


The  Chapters  to  Follow 


Chapter  II  briefly  examines  a  single  sampling  plan  for  variables 
inspection  and  then  presents  the  development  of  Wald's  SPR  test  for 
variables  inspection.  Also  discussed  are  the  methods  used  to 
calculate  the  OC  curve  and  average  sample  number  (ASN)  curve  of  the 
SPR  plan.  Finally  various  truncation  rules  and  their  effect  on  the 
plan's  discriminatory  power  are  presented.  The  SNPR  procedure  and  a 
heuristic  method  of  estimating  the  procedures  arguments  are  presented 
in  Chapter  III  followed  by  a  comparison  of  SPR  and  SNPR  procedures  in 
Chapter  IV. 


CHAPTER  II 


HISTORICAL  PERSPECTIVE 


Single  Sampling  Plan 


A  single  sampling  plan  for  variables  inspection  is  relatively 
simple  to  design.  Given  X1^,  ^'lTPD*  a’  6  and  0,>  a  sample  size  N 


and  an  acceptance  limit  X,  are  determined  to  establish  the  plan.  The 

a 


derivation  of  the  plan  is  as  follows.  If  a  lot  actually  has  a  mean 
grade  of  AQL  quality,  means  of  samples  of  N  from  the  lot  will  have 


a  normal  distribution  with  mean  X ' Aqj_  and  a  standard  deviation  of 


a'//N.  Therefore,  if  a  is  the  probability  of  rejecting  a  lot  with 


mean  X 1 ,  then 


y  _  Y 1 

*a  A  AQL  _ 


7/N 


(1) 


Similarly,  a  lot  having  a  mean  grade  of  X'^-j-pQ  and  a  probability 
acceptance  B  requires 


x  -  x' 
a  A  LTPD 


'//N 


=  -Z, 


(2) 


When  X'AqL,  X'^-j-pg,  a,  B  and  a'  are  defined,  we  have  two  equations  in 
Xa  and  N.  The  simultaneous  solution  of  which  yields  values  of  X,  and 


The  SPR  plan  for  variables  [  1  ]  is  one  which  requires  the 
sequential  observation  of  a  variable  X  which  is  normally  distributed 
with  mean  X'  and  standard  deviation  a'.  The  decision  variable  is 


defined  as: 


T(n)  =  l  X. 
i=l 


(3) 


where  n  is  the  number  of  observations.  The  operation  of  the  plan  is 
described  by  analyzing  the  sequence  {T(n)h  Consider  the  point 
T(1 )  =  X-j  where  X-j  is  the  first  sample  observation.  If  T(l)  falls 
on  or  above  an  upper  bound  Tu(l)  and  low  values  of  X  are  desirable, 
the  lot  is  rejected.  If  T(l)  falls  on  or  below  a  lower  bound  T^(l), 
the  lot  is  accepted.  If,  however,  T£(l)  <  T(l)  <  Tu(l),  another 
observation  is  taken.  Then  T( 2)  =  Xi  +  X2  is  evaluated  in  a  similar 
manner,  and  so  on  until  a  final  course  of  action  can  be  determined. 

The  decision  variable  T(n)  used  for  testing  the  mean  of  a  normal 
2 

distribution  with  o'  known  is  derived  from  Wald's  sequential 


probability  ratio  test  [9].  The  test  is  given  as  follows.  Two 
points  (9g,  a)  and  (e^,  8)  are  selected  where  a  is  the  probability  of 
rejecting  X ‘  =0 q  and  is  the  probability  of  accepting  X '  =0 -j .  The  SPR  of 
the  random  sample  (X,,  X2>  ...»  X„)  from  the  normal  distribution  is 


defined  as 


!  1  l  (X,  -  e  )2 

1  e - i=l  1  1 

(2ir)n/2a'n  2°‘2 

SPR  =  - 

_j _  _  _L  X  (Xi "  e°)2 

(2Tf)n/2CT,n  2°'2 


(4) 


In  theory,  the  ratio  is  computed  at  each  stage  of  the  inspection 
and  additional  observations  are  taken  as  long  as 

B  <  SPR  <  A  (5) 

where  B  and  A  are  appropriately  defined  constants.  Inspection  is 

terminated  with  acceptance  of  the  lot  if 

SPR  <  B  (6) 

Inspection  is  terminated  with  rejection  of  the  lot  if 

SPR  >  A  (7) 


Selecting  A  =  (l-e)/a,  B  =  B/(l-a),  taking  logarithms  of  the  inequal 
ties  and  simplifying,  (5),  (6)  and  (7)  become 


T(n)  s  X  xi  r^e-  log  +  n 

i=l  1  el  e0  1  °  < 


1-3  ,  „  90+el 


T(n)  -  l  X.  >  log  ^  +  n  ^ 


respectively.  The  decision  variable  therefore  becomes  T(n)  rather 
than  the  SPR. 

By  using  X ' ^qj_  -  eQ  and  ^'LTPD  =  el»  the  for  the  Plan  from 

(9)  and  (10)  become 


Tu(n)  =  n1  +  S  n 


T^n)  =  hQ  +  S  n 


where  the  common  slope  formed  by  the  limit  lines  is  given  by 


Y  l  +  y  1 

A  AQL  A  LTPD 


and  the  intercept  of  the  lines  when  n=0  is  given  by 


.  -  g  log  B 

0  "  S'  S' 

A  LTPD'  AQL 


h  =  _ 1 Q3L  A 

1  Y 1  _  Y ' 

A  LTPD  A  AQL 


(15) 


The  SPR  PC  and  ASN  Curves 


Wald  19]  shows  that  on  an  OC  curve  for  X'  =  -®,  ^ ' AQL ’  ^'lTPI 

^'aQL  +  *LTPD^2'  and  +  the  va1ues  of  L(X')  (the  pa  of  a  lot 
whose  mean  equals  X')  are  as  follows: 


L(-)  »  1  L(X'aql)  =  1-a 


log  A 

log  A  -  log  B 


(16) 


UX'ltrd)  =  6  L(-)  =  0 

An  approximate  general  equation  is  given  by 
Ah  -  1 

L(X')*\ - \  (17) 

An  -  Bn 


where 


h  = 


X 1  +  Y1  -  ?X' 

A  AQL  A  LTPD  ^ 

V  -  x' 

A  LTPD  A  AQL 


(18) 


Wald  also  derived  an  approximate  formula  for  the  ASN  curve  which  is 
given  as 

_  ..  h.  +  L(r)  (h„-h.) 


Using  the  values  defined  in  Chapter  I,  the  number  of  observations 
and  acceptance  and  rejection  limits  are  tabulated  in  Table  1  along 
with  a  sample  set  of  observations.  The  same  test  plan  is  graphically 
illustrated  in  Figure  3.  The  corresponding  OC  and  ASN  curves  are 
shown  by  Figure  4  and  Figure  5  respectively.  As  indicated,  the  plan 
would  terminate  with  acceptance  of  the  lot  at  n  =  8. 

Truncation  of  the  SPR  Plan 

As  mentioned  previously,  the  main  disadvantage  of  the  SPR  plan  is 
that  there  is  no  effective  upper  limit  to  the  number  of  items  that  are 
required  for  inspection  before  a  decision  to  accept  or  reject  a  lot 
can  be  made.  The  sample  size  of  an  SPR  plan  is  a  random  variable  whose 
value  will  occasionally  be  quite  large.  Truncation  of  the  procedure 
is  used  to  force  a  decision  prior  to  or  at  n  =  n'  where  n'  is  an  a- 
priori  fixed  integer.  The  following  is  a  common  rule  for  accepting 
or  rejecting  a  particular  lot  if  a  decision  has  not  been  reached  for 
n  <  n'  with  the  regular  sequential  procedure.  If  T(n')  > 

(T£{n‘)  +Ta(n'))/2  the  lot  is  rejected,  and  if  T(n')  <  (T£(n')+Tu(n'))/2 
the  lot  is  accepted. 

Several  rules  of  thumb  have  been  suggested  for  establishing  an 
effective  n'.  Guild  [5]  uses  n'  equal  to  twice  the  maximum  ASN  at 
X ' aql  or  *'ltpd*  Grant  and  Leavenworth  [4]  set  n1  equal  to  three  times 
N;  the  number  of  samples  required  by  a  single  sampling  plan  with  the 


Table  1.  Acceptance  and  Rejection  Limits  of  Sample 
Problem  SPR  Plan  With  Sample  Observations 


CUMULATIVE  SAMPLE  SIZE,  n 


Figure  3.  Plot  of  Sample  Problem  SPR  Plan  with  Observations 


same  level  of  protection.  Wald  [9  ]  indicates  that  with  n'  put  as 
high  as  three  times  the  maximum  expected  value  as  indicated  by  the 
ASN  curve,  the  effect  of  truncation  on  the  OC  curve  is  negligibly 
small  since  the  probability  is  nearly  one  that  the  regular  sequential 
procedure  will  terminate  for  n  <  n'.  For  the  Wald  method  of  trunca¬ 
tion,  the  value  of  n'  to  be  used  in  the  sample  problem  is  calculated 
to  be  30.  Compared  to  13,  the  number  of  samples  required  by  the 
single  sampling  plan  to  yield  the  same  level  of  protection,  30  seems 
relatively  large  and  possibly  prohibitively  expensive  if  the  cost  of 
inspection  is  high. 

Table  2  is  extracted  from  Wald's  "Sequential  Analysis"  [9]  and 
illustrates  the  effect  of  truncation  on  the  plan  for  different  values 
of  a  and  6.  If  the  plan  is  based  on  the  values  of  a  and  6  shown,  but 
a  decision  is  made  at  n1  even  when  the  regular  sequential  plan  requires 
a  continuation  of  the  process,  the  realized  values  a(n')  and  B ( n * ) 
will  not  exceed  the  tabular  entries.  The  table  relates  to  the  test  of 
the  mean  of  a  normally  distributed  variate,  the  difference  between  the 
null  and  alternate  hypothesis  adjusted  for  each  pair  (a,  b)  so  that  the 
number  of  trials  required  by  the  single  sampling  plan  of  strength 
(a,  b)  is  1000.  Notice  that  as  a  and  6  are  relaxed,  the  values  of  a(n') 
and  B(n'),  even  when  n'  is  increased  to  twice  N,  vary  by  15%  when 
a  =  .05  and  B  =  .05.  Wald's  opinion  is  that  the  upper  limits  given  in 
the  table  are  considerably  above  the  true  values  of  a(n')  and  B(n') 
when  n'  is  not  much  higher  than  N.  However,  any  difference  between  the 
pairs  (a ,  b)  and  (a(n'),  s (n ' ) )  may  be  intolerable.  A  more  reasonable 


Table  2.  Effect  on  Risks  of  Error  of  Truncating  a 
Sequential  Analysis  at  a  Predetermined 
Number  of  Trials 


a  =  .01  and  8  =  .01 

a  =  .01  and  8  =  .05 

a  =  .05  and  8  =  .05 

Number 

of 

Trials 

Upper 
bound  of 
effective 

Upper 
bound  of 
effective 

Upper 
bound  of 
effective 

Upper 
bound  of 
effective 

Upper 
bound  of 
effective 

Upper 
bound  of 
effective 

a(n' ) 

B  (n 1 ) 

a(n') 

B(n') 

a(n' ) 

B(n') 

1000 

.020 

.020 

.033 

.070 

.095 

.095 

1200 

.015 

.015 

.024 

.063 

.082 

.082 

1400 

.013 

.013 

.019 

.058 

.072 

.072 

1600 

.012 

.011 

.016 

.055 

.066 

.066 

1800 

.011 

.010 

.014 

.053 

.062 

.062 

2000 

.010 

.010 

.012 

.052 

.058 

.058 

2200 

.010 

.010 

.012 

.051 

.056 

.056 

2400 

.010 

.010 

.011 

.051 

.055 

.055 

2600 

.010 

.010 

.011 

.051 

.053 

.053 

2800 

.010 

.010 

.010 

.050 

.053 

.053 

3000 

.010 

.010 

.010 

.050 

.052 

.052 
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plan  would  be  one  which  guarantees  the  same  level  of  protection  with 
a  commensurate  savings  in  required  sample  size. 


CHAPTER  III 


SEQUENTIAL  NON-PROBABILITY  RATIO  SAMPLING  PLAN 

The  SNPR  plan  is  an  acceptance  sampling  procedure  bounded  by 
two  lines  which  intersect  at  n  =  n',  where  n'  is  an  a-priori  fixed 
integer.  Although  the  procedures  for  SPR  and  SNPR  plans  are  similar, 
the  SNPR  plan  requires  no  truncation  rules.  As  mentioned  in  the 
previous  chapter,  truncation  of  the  SPR  plan  makes  the  OC  curve  of 
the  plan  less  discriminating.  Therefore,  the  implicit  property  of 
the  SNPR  plan  may  be  desirable  if  it  can  provide  the  same  level  of 
protection  as  the  non-truncated  SPR  plan. 

The  General  Approach 

For  the  SNPR  plan,  T£(n)  and  ^(n)  are  constructed  which  converge 
at  some  maximum  number  of  units,  n'.  The  decision  variable  for  the 
plan  is  the  same  as  that  defined  in  equation  (1)  for  the  SPR  plan. 
Similarly,  the  operation  of  the  plan  is  the  same.  If  no  decision  is 
made  prior  to  sampling  the  n'th  unit,  the  lot  is  rejected  if 
T(n')  >  n ' X ‘  and  accepted  if  T(n')  <  n'X1.  The  plan  is  graphically 
illustrated  in  Figure  6. 

T^  (1 )  and  T  ( 1 )  are  determined  from 


CUMULATIVE  SAMPLE  SIZE 

Figure  6.  Plot  of  an  SNPR  Plan 
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so  that  the  probability  of  acceptance  (P  1  and  the  probability  of 

a 

rejection  (P  )  are  both  equal  to  a'  for  X'  at  n  =  1.  At  n  =  n' 

Tu(n')  =  Tjn')  =  n'X'  (22) 

The  convergence  by  the  limits  on  n'X1  indicates  a  step  by  step  reduc¬ 
tion,  d,  in  the  size  of  the  interval  Z  ,a'.  This  is  shown  by 

a 

d  =  Za , a 7 ( n 1  - 1 ) ,  n*  >  0  (23) 

Therefore,  the  calculation  for  the  upper  limit  becomes 

Tu(n)  =  nX'  +  Zo,o'  -  (n-1 )  Zo,<y'/(n'-l) 
or  (24) 

T^(n)  =  nX'  +  (n'-n)Za,o7(n'-l),  1  <  n  <  n* 

Similarly, 

T  (n)  =  nX'  -  (n'-n)Z  ,o7(n'-l)t  1  <  n  <  n'  (25) 

X»  Cl 

X'  and  Z  ,  for  the  SNPR  plan  must  now  be  established  so  that  limits 
a 

which  provide  the  same  level  of  protection  as  the  SPR  plan  may  be 
calculated. 


Sample  Illustration 


Using  the  example,  if  T  (1)  and  T^(l)  are  arbitrarily  set  to  the 
same  values  as  in  the  SPR  plan,  X'  and  ZqI  from  equations  (20)  and  (21) 


become 


Tu(1)  +  T£  1 ^  _  316.55  +  131.46 


a 1  = 


=  3.0840,  a'  =  .0010 


Table  3  shows  the  proposed  plan  for  n'  =  10  (10  is  the  maximum  ASN  of 

the  SPR  plan)  while  Figure  7  shows  the  plot. 

a  and  8  for  the  SNPR  plan  are  calculated  in  the  following  manner, 
n  o  o 

Since  T(n)  =  l  X.  and  X  is  N(X',  a'^),  T(n)  is  N(nX',  no'^)  and 
i=l  1 

a(T(n))  =  /no'.  Therefore,  the  normal  deviate 


Z  =  (Tu(n)  -  nX' )//na' 


can  be  used  to  solve  for  a.  Table  4  lists  the  results  of  calculations. 
Using  the  emotion 


a  =  a-|  +  (1 -a-| -8-|  +  (1 -a.j -B-| )  ( l-ag-Bg)0^  +  •••  + 

...  (l-cxg-ag)a^g 

where  aR  =  P(Z>Za  | X‘  =X‘ AQL)  and  Bk  =  P(Z<Zg  | X*  =X* LTp[)) ,  a  is  cal- 
n  n 

culated  to  be  .005.  Similarly  8  is  calculated  to  be  .450.  Figure  8 
shows  the  0C  curves  for  the  corresponding  SPR  and  SNPR  plans. 

The  arbitrary  values  of  X'  and  Z^,  chosen  for  the  SNPR  plan 
yield  an  0C  curve  which  is  obviously  less  discriminating  than  that 
of  the  SPR  plan.  However,  altering  the  values  of  X',  Z^,  and  n' 
changes  the  values  of  a  and  8.  It  will  be  shown  that  properly  chosen 
values  of  X',  Za,  and  n'  will  yield  an  0C  curve  for  the  SNPR  plan 
which  discriminates  as  well  as  the  0C  curve  of  the  SPR  plan. 


Table  3.  Acceptance  and  Rejection  Limits  of  Sample 
Problem  SNPR  Plan  (X1  =  224.01,  Z^,  =  3.0 
and  n'  =  10)  With  Sample  Observations 


X(n) 

T(n) 

V"> 

Tu(n) 

201 

201 

131.49 

316.53 

206 

407 

365.78 

530.26 

199 

606 

600.07 

743.99 

204 

810 

834.36 

957.72 

1068.65 

1171.45 

1302.94 

1385.18 

1537.23 

1598.91 

1771.52 

1812.64 

2005.81 

2026.37 

2240.10 

2240.10 

Decision 


SAMPLE 

ACCEPT 


SUM  OF  OBSERVATIONS,  T(n) 


Table  4.  Calculations  Used  to  Solve  for  ot  and  3  in  Sample  Problem  SNPR  Plan 
(X *  =  224.01  ,  Z  ,  =  3.0840  and  n1  =  10 


.0000  +  (.9887)1.0010)  ♦  (.9887)(.7801){.0027)  +  ...  +  ...  (0)  .0056 


PROBABILITY  OF  ACCEPTANCE 
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LOT  x' 


Figure  8.  OC  Curves:  SNPR  Plan  (Sample  Problem;  X1  =  224.01, 

Za,  =  3.0840  and  n'  -  10)  and  SPR  Plan  (Sample  Problem) 


Preliminary  Analysis 

Holding  n'  and  ZqI  constant,  X1  is  altered  to  examine  the  effect 
on  a  and  8.  Table  5  lists  the  results.  It  is  observed  that  as  X' 
decreases,  a  increases  and  8  decreases.  This  is  intuitively  reasonable 
because  altering  X'  relocates  the  envelope  formed  by  the  acceptance  and 
rejection  limits.  A  downward  relocation  of  the  envelope,  caused  by 
decreasing  X',  increases  the  rejection  region.  If  the  rejection  region 
is  increased,  the  probability  of  rejecting  acceptable  material  is  in¬ 
creased,  therefore,  a  is  increased.  Similarly  8  decreases  as  X' 
decreases.  This  is  graphically  illustrated  in  Figure  9.  OC  curves  for 
the  same  plans  are  compared  to  the  SPR  plan's  OC  curve  in  Figure  10. 

A  strong  similarity  exists  between  the  OC  curves  of  the  SPR  plan  and 
the  SNPR  plan  with  X'  =  214.05.  This  particular  value  of  X’has  not  been 
chosen  arbitrarily.  It  is  the  value  of  X  .  the  acceptance  limit  for  the 

u 

single  sampling  plan  discussed  in  Chapter  II.  X  will  be  used  in  the 

cl 

development  of  a  heuristic  procedure  which  establishes  X'  for  an  SNPR 
plan  equivalent  to  the  SPR  plan. 

Varying  Z^,  and  holding  X'  and  n'  constant  effects  the  values  of 
a  and  8  as  shown  in  Table  6.  As  Z^,  increases,  a  and  8  both  decrease. 
Again  this  is  intuitively  reasonable  because  increasing  ZqI  widens 
the  envelope  formed  by  the  acceptance  and  rejection  limits.  The 
widening  decreases  both  the  acceptance  and  rejection  regions,  as  shown 
in  Figure  11,  thereby  decreasing  both  a  and  8.  Figure  12  shows  the 


corresponding  OC  curves.  Again  a  similarity  exists  between  the  OC 


SUM  OF  OBSERVATIONS,  T(n) 


Figure  9,  Plot  of  Sample  Problem  SNPR  Plans  (I  ,  =  3.0840, 
n'  =  10  and  X'  =  214.05,  224.01  and  230.00) 


OBSERVAT 


CUMULATIVE  SAMPLE  SIZE,n 


Figure  11.  Plot  of  Sample  Problem  SNPR  Plans  (X1  =  214.05, 
n'  =  10  and  Z  ,  =  2.0000,  3.0840  and  4.0000) 


TY  OF  ACCEPTANC 


Table  6.  Effect  on  a  and  6  of-v^lte^^n,9nZ,a,  jSf?o) 
Problem  SNPR  Plan:  X*  =  214.05,  n  iu; 


z  . 

a 

Tt(1) 

V” 

T(n') 

2.000 

154.05 

! 

274.04 

2140.50 

3.084 

121.53 

306.57  ! 

2140.50 

4.000 

94.05 

334.04 

2140.50 
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OC  Curves:  SNPR  Plans  (Sample  Problem; 
l'  =  214.05,  n'  =  10  and  7^,  -  2.0000, 
3.0840  and  4.0000)  and  SPR  Plan 


curves  of  the  SPR  plan  and  the  SNPR  plan  with  X'=  214.05  and  Z^,  =  3.0840 
These  values  are  used  in  the  development  of  a  heuristic  which  for  a 
specified  value  of  n'  establishes  X'  and  Z^,  for  an  SNPR  plan  which 
provides  the  same  level  of  protection  as  the  SPR  plan. 


Development  of  the  Heuristic 

Three  important  facts  mentioned  in  the  previous  section  are  the 
cornerstones  of  the  heuristic  procedure: 

(1)  The  OC  curve  of  an  SNPR  plan  with  X1  equal  to  Xa  of  the  single 

a 

sampling  plan,  Z^,  equal  to  the  Zq,  obtained  from  the  limits 
at  n=l  of  the  SPR  plan  and  n'  equal  to  the  maximum  ASN  of  the 
SPR  plan  closely  approximates  the  OC  curve  of  the  SPR  plan. 

(2)  Increasing  the  value  of  X'  in  an  SNPR  plan  decreases  the  value 
of  a  and  increases  8.  Decreasing  X'  has  the  opposite  effect. 

(3)  Increasing  the  value  of  Z^,  in  an  SNPR  plan  decreases  both  a 
and  6.  Decreasing  Za,  has  the  opposite  effect. 

Using  this  information,  Z^,  and  X'  are  altered  while  holding  n' 
constant  at  10.  By  setting  X'  =  214.24  and  Z^,  =  3.0100,  the  calculated 
values  of  a  and  e  are  .05  and  .10  respectively.  If  n  is  reduced  to  9 
(a  10%  sampling  reduction),  the  values  of  X'  and  Z^,  need  to  yield 
a  =  .05  and  8  =  .10  are  214.22  and  3.2400  respectively.  Table  7  lists 
the  necessary  values  of  X1  and  Z ,  for  different  values  of  n'.  Figure  13 
shows  the  plots  of  two  plans  along  with  a  plot  of  the  SPR  plan.  Notice 
that  as  n'  decreases,  the  plan  compensates  by  being  less  discriminating 
in  its  earlier  stages.  Computational  experience  indicates  that  SNPR 


Table  7.  Values  of  X'  and  Za,  Necessary  to  Obtain 
a  =  .05  and  B  =  .10  for  Varying  Values  of 
n'  in  the  Sample  Problem 


CUMULATIVE  SAMPLE  SIZE  ,  n 


Figure  13.  Plots  of  Sample  Problem  SNPR  Plans  (n1  =  7  and 
n'  =  10)  and  Sample  Problem  SPR  Plan 


plans  with  n'  less  than  70%  of  the  SPR  plan's  maximum  ASN  cannot  yield 

an  OC  curve  which  approximates  that  of  the  SPR  plan,  regardless  of  the 

values  of  X'  and  Z  ,  chosen. 

a 

Using  the  estimates  of  X1  and  Z^,  described  at  the  beginning  of 
this  section,  a  more  exacting  OC  curve  can  be  obtained  using  the 
following  steps: 

(1)  If  n1  is  chosen  other  than  the  maximum  ASN  of  the  SPR  plan 
(Nmax),  a  modified  estimate  of  Z^,  is  used.  This  value  is 
obtained  by  increasing  (decreasing)  the  initial  estimate  by 
the  percentage  decrease  (increase)  in  Nmax: 

Za,(n')  =  Z^, (Nmax)(l+(l-n'/Nmax) )  (26) 

(2)  Examine  the  values  of  a  and  8  obtained  by  using  the  estimated 
parameters  in  the  SNPR  plan.  Determine  if  they  must  be 
increased  or  decreased  to  match  the  objective  values  of  a  and 
6  from  the  SPR  plan. 

(3)  If  both  a  and  8  must  be  increased  (decreased)  to  meet  the 
objective  values  of  a  and  8  from  the  SPR  plan,  decrease 
(increase)  Z  ,. 

a 

(4)  If  one  (a  or  b)  must  be  increased  and  the  other  decreased  to 
meet  the  objective  values  of  a  and  8  from  the  SPR  plan, 
alter  X1  accordingly. 

(5)  Repeat  steps  2  through  4  until  the  values  of  a  and  8  from 
the  SNPR  test  are  within  the  desired  limits. 

Computational  experience  with  the  heuristic,  as  coded  in  Appendix  A, 
reveals  that  usually  only  two  or  three  iterations  of  the  procedure  are 
needed  to  obtain  a  and  8  within  a  tolerance  of  .001.  Appendix  B  contains 


several  SNPR  plans  which  were  generated  using  the  coded  heuristic. 

The  associated  problem  parameters  were  varied  through  specific  ranges 
to  demonstrate  the  code's  flexibility.  In  each  of  the  twelve  plans, 
the  exact  objective  values  of  a  and  6  were  met. 

Using  the  sample  problem,  the  following  is  an  example  of  the 
heuristic  procedure  for  n'  =8. 

(1)  Initial  Estimates 

(Step  1)  V  =  214.05  Z  ,  =  3.0849 

Z  ,(8)  =  3.0849  (1+0-8/10))  =  3.7019 

(2)  With  X'  =  214.05,  Z  ,  =  3.7019  and  n'  =  8 

(Step  2)  a  =  .0507  8  =  .0965 

(3)  Since  a  is  above  the  objective  value  of  .05  and  8  is  below 
(Step  4)  the  objective  value  of  .10,  X'  is  increased.  (The  rule  for 

altering  X'  is  to  first  alter  it  by  o'/lOO.  After  recalcu¬ 
lating  the  values  of  a  and  8,  determine  the  proportionate 

adjustment  to  X1  which  yields  a  and  8  closest  to  the  objec- 

,  * 
tive  values. ) 

With  X'  =  214.35,  Z  ,  =  3.7019  and  n'  =  8 

a 

a  =  .0474  6  =  .1023 

a  proportionate  adjustment  yields  X'  =  214 
With  X'  =  214,  Z  ,  =  3.7019  and  n'  =  8 
a  =  .0495  8  =  .0993 

(4)  With  a  and  8  within  .001  of  the  objective  values  of  .05  and 
.10,  respectively,  the  heuristic  procedure  is  terminated. 

ic 

The  rule  for  altering  Z  ,  is  to  increase  (decrease)  the  value  of  Z  , 
accordingly,  recalculate  the  values  of  a  and  8,  then  determine  the 
proportionate  adjustment  to  Z  ,  which  yields  a  and  8  closest  to  the 
objective  values.  a 


CHAPTER  IV 


COMPARISON  OF  SPR  AND  SNPR  PLANS 

Because  it  reveals  the  discriminating  power  of  a  sampling  plan, 
the  OC  curve  is  useful  in  evaluating  alternate  sampling  plans.  In  the 
previous  chapter  a  method  to  closely  match  the  OC  curves  of  the  subject 
plans  was  presented  and  therefore  voids  any  comparisons  based  on  that 
feature.  However,  situations  may  exist  where  it  is  advantageous  to 
quickly  establish  a  sequential  sampling  plan  which  guarantees  truncation 
after  relatively  few  samples  have  been  taken,  even  if  the  cost  of  trun¬ 
cation  is  a  slightly  altered  OC  curve.  Two  methods  are  available  to 
establish  a  "fast  truncated"  sequential  sampling  plan.  They  are  the 
truncated  SPR  plan  described  in  Chapter  II  and  the  SNPR  plan  which  uses 
the  initial  estimates  of  X'  and  Z  ,  (5L  from  the  single  sampling  plan 

Cl  a 

and  Z^,  derived  using  the  initial  acceptance  and  rejection  limits  of 
the  SPR  plan)  discussed  in  Chapter  III. 

Four  sampling  plans  are  presented  in  Table  8.  The  plans  are  the 
truncated  SPR  plan  (n '  =  10),  the  non-truncated  SPR  plan,  the  SNPR  plan 
(n1  =  10)  generated  using  the  coded  heuristic  procedure  in  Appendix  A 
and  the  SNPR  plan  (n1  =  10)  generated  using  the  initial  es  imates 
mentioned  above.  The  first  plan  (truncated  SPR)  has  a(10)  =  .119  and 
B ( 1 0 )  =  .156  (these  values  are  upper  limits  for  a(n')  and  B ( n ' )  com¬ 
puted  via  equations  developed  by  Wald  [9]).  While  this  is  a  "fast 
truncated"  sequential  plan,  the  resulting  distortion  of  the  OC  curve 
leaves  the  practicality  of  this  plan  highly  suspect.  The  fourth  plan 
(SNPR  generated  using  initial  estimates)  has  a  =  .050  and  8  =  .095. 
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The  0C  curve  of  this  plan,  although  not  exactly  matched  to  the  0C  curve 
of  the  non-truncated  SPR  plan,  may  reveal  sufficient  protection  that  it 
would  not  be  necessary  to  perform  the  heuristic  procedure  to  obtain  a 
more  exacting  0C  curve. 

If  exact  discriminating  power  of  a  plan  is  an  absolute  requirement, 
Duncan  [1]  suggests  the  relative  efficiency  of  plans  with  matching  0C 
curves  may  be  determined  by  examining  the  amount  of  inspection  required 
by  each  plan.  The  ASN  curve  of  the  SPR  plan  shown  in  Figure  5  is 
characteristic  of  all  SPR  plans.  It  can  be  seen  that  decisions  (generally 
to  accept)  are  quickly  made,  on  the  average,  when  good  quality  lots  are 
submitted  and  that  decisions  (generally  to  reject)  are  quickly  made,  on 
the  average,  when  bad  lots  are  submitted.  Lots  of  mediocre  quality 
require  more  inspection,  on  the  average  [7],  The  phrase,  "on  the 
average"  is  stressed  because  the  required  sample  size  is,  as  mentioned 
previously,  a  random  variable  capable  of  taking  on  values  well  above 
the  average. 

Although  no  analytic  function  has  been  developed  which  describes 
or  approximates  an  SNPR  plan's  ASN  curve,  it  is  relatively  easy  to 
envision  its  rudimentary  shape.  It  is  known  that  values  on  the  curve 
cannot  exceed  n'  but  will  most  likely  take  on  values  well  below  n'. 
Therefore,  n1  represents  an  upper  limit  for  the  curve.  In  addition. 


Figure  13  are  plots  of  three  sampling  plans  developed  from  the  sample 
problem.  In  the  initial  stages,  n  =  1  through  4,  the  SNPR  plan 
( n '  =  10)  is  as  discriminating  as  the  SPR  plan  while  the  SNPR  plan 
(N‘  =  7)  is  less  discriminating.  This  suggests  that  the  ASN  curve  of 
the  SNPR  plan  (n1  =  10)  closely  resembles  that  of  the  SPR  plan  at  the 
tails  of  the  ASN  curve.  Conversely,  the  ASN  curve  of  the  SNPR  plan 
(n1  =  7)  is  above  that  of  the  SPR  plan  at  the  tails  of  the  ASN  curve. 

If  it  is  assumed  that  incoming  lots  of  mediocre  quality  force  the 
SNPR  plan's  ASN  curve  to  approach  n',  the  ASN  curves  in  Figure  14  are 
rough  approximations  of  what  can  reasonably  be  expected. 

Examining  the  curves  in  Figure  14,  an  expected  feature  of  the  SNPR 
plan  becomes  evident.  As  n1  decreases,  so  does  the  maximum  average 
sample  number.  However,  the  curve  becomes  less  peaked;  the  average 
sample  number  increases  at  values  of  X1  which  indicate  lots  of  very 
good  or  very  bad  quality.  At  this  point,  a  trade-off  must  be  made  when 
selecting  which  plan  to  use.  If  one  is  reasonably  sure  that  the  values 
of  X'  are  clustered  near  one  value,  the  selection  of  a  plan  can  be 
made  by  comparing  ASN  curves  and  selecting  the  plan  with  the  lowest 
ASN  at  that  point.  If  the  quality  of  incoming  lots  is  unknown,  the 
cost  of  inspection  versus  the  number  of  inspections  trade-off  analysis 


AVERAGE  SAMPLE  NUMB 


LOT  x' 


Figure  14.  ASN  Curves:  S'!PR  Plans  (Sample  Problem  n1  =  7,  10) 
and  SPR  Plan  (Sample  Problem) 
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CHAPTER  V 


SUMMARY 


An  SPR  plan  can  materially  reduce  the  amount  of  required  inspection. 
Studies  have  shown  that  the  average  decrease  in  sample  size  is  often 
near  50%  when  compared  with  the  sample  size  of  a  comparable  single 
sampling  plan  [4  1.  An  SNPR  plan  can  achieve  and  better  the  sampling 
reduction  obtained  through  the  use  of  an  SPR  plan.  The  SNPR  plan 
requires  no  truncation  rules  and  maintains  the  integrity  of  its 
associated  OC  curve  while  truncation  of  the  SPR  plan  makes  the  plan 
less  discriminating. 

Several  disadvantages  of  sequential  plans  have  been  noted.  The 
most  severe  are  the  cost  of  administering  the  plan  and  the  lack  of 
availability  of  information  concerning  prevailing  levels  of  quality  in 
each  lot.  The  administrative  cost  of  an  SNPR  plan  can  be  substantially 
reduced  through  the  use  of  the  computer  code  of  the  heuristic  algorithm 
developed  in  Chapter  III.  No  bulky  tables  and  complicated  charts  need 
be  used,  only  a  terminal  display.  Acceptance  and  rejection  limits  would 
be  displayed  on  a  screen  in  conjunction  with  input  inspection  variables. 
The  savings  in  inspection  costs  and  time  could  possibly  offset  any  loss 
due  to  unavailability  of  lot  quality  data.  The  ultimate  extension 
would  be  the  use  of  the  sampling  plan  in  an  automated  environment  where 
inspection  is  performed  by  machines. 

Whatever  its  use,  the  "NPR  plan  for  variables  inspection  is  a 
flexible  plan  which  in  many  instances  outperforms  other  sampling  plans. 
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Its  applications  are  not  limited  to  one-sided  inspection  and  further 
research  should  be  able  to  produce  modifications  which  enable  it  to  be 
expanded  to  two-sided  inspection  (inspection  with  both  upper  and  lower 
quality  limits).  It  is  hoped  that  as  a  result  of  this  research, 
serious  consideration  will  be  given  to  this  plan  so  as  to  reap  the 
potential  savings  that  can  be  a  by  product  of  its  use. 


APPENDIX  A 


SEQUENTIAL  NON-PROBABILITY  RATIO  (SNPR)  AND 
SEQUENTIAL  PROBABILITY  RATIO  (SPR) 

TEST  PLAN  GENERATOR 

The  following  program  provides  a  SNPR  test  plan,  a  SPR  test  plan 
or  both  for  a  given  set  of  input  parameters.  It  also  calculates  a 
specified  number  of  points  for  the  associated  OC  curve(s)  of  the 
requested  plan(s).  Additionally  a  specified  number  of  points  are 
calculated  for  the  ASN  curve  of  the  SPR  plan  when  that  plan  alone  is 
requested. 

The  program  has  three  formats: 

(1)  the  SPR  test  plan; 

(2)  the  SNPR  test  plan; 

(3)  the  SPR  and  SNPR  test  plans. 

The  following  is  a  brief  description  of  each  format's  input  parameters 
and  output.  The  input  format  is  "format  free,"  requiring  only  commas 
or  spaces  between  input  parameters  for  each  line  of  data.  Before  input 
data  can  be  read  for  any  particular  format,  a  numerical  format  code 
must  be  specified  on  the  first  data  card.  The  format  specification 
codes  are: 

CODE  FORMAT 

(1)  SPR  Test  Plan 

(2)  SNPR  Test  Plan 

(3)  SPR  and  SNPR  Test  Plans 


CARD ( 1 ) :  INPUT  FORMAT  CODE 
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0.6299 

4 

9.22 

in 

* 

2653.95 

4 

3079. 05 

212.71 

4 

0.55  16 

4 

9.  3e 

IS 

• 

3106.95 

* 

3291.55 

219.17 

4 

0.9779 

4 

9.  91 

Ifi 

4 

33 ie.ss 

4 

3509.05 

215.62 

4 

0.9051 

4 

9.  30 

17 

4 

3531.95 

4 

3716.55 

217.08 

4 

C . 33  75 

4 

9.07 

IP 

4 

3793.55 

4 

3929.05 

218.59 

4 

0.2766 

4 

8.  75 

19 

4 

3956.95 

4 

9191.55 

220.00 

4 

0.2235 

* 

8.36 

20 

4 

9168.95 

4 

9359.05 

221. 96 

4 

C.  1786 

4 

7.93 

4  4 

222 .92 

4 

0.1919 

4 

7.  98 

229. 37 

4 

0.1111 

4 

7.  03 

225. 83 

4 

o.oeee 

4 

6.60 

227.29 

4 

C . 0675 

4 

6.  19 

229.75 

4 

0.0529 

* 

5.81 

230  .  2  1 

4 

1.0905 

4 

5.  95 

231  .67 

4 

0.0313 

4 

5.  13 

233.  1  2 

4 

0.0292 

4 

9.  P3 

234.58 

4 

C  .0  186 

4 

9.  66 

* 

44 

236 . OU 

*4*4*444 

* 

* 

0.0193 

4 

*  * 

9.  31 

Figure  15.  Program  Output  Listing:  Format 
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LOTXBAR  --  values  of  X1  on  x  axis  of  CO  and  ASN 
PROB  ACC  --  probability  of  acceptance  for  X' 

ASN  --  average  sample  number  of  X' 

The  following  are  the  data  cards  necessary  for  format  1  input: 

CARD  (2):  X'AQL.  X'LTp[),  o' 

CARD  (3):  ct,  B 

CARD  (4):  NTRUNK,  NTERM 

Format  2:  SNPR  Test  Plan 

The  input  format  for  the  SNPR  test  plan  requires  specification  of: 


(1) 

x 1 

X  AQL 

(2) 

x 1 

A  LTPD 

(3) 

a' 

(4) 

n' 

(5) 

a 

(6) 

0 

(7) 

MAX 

(8) 

NTERM 

All  but  two  of  the  above  terms  have  been  described  previously,  n'  is 
the  maximum  number  of  items  to  be  sampled.  MAX  is  the  maximum  number 
of  iterations  to  be  performed  by  the  heuristic  algorithm.  This  value 
should  initially  be  set  to  5  and  adjusted  depending  on  results  obtained. 

Figure  16  shows  an  example  of  an  input  list  and  the  associated 
output.  In  the  output  listing,  a  header  echoes  the  input  and  lists 
ZA  and  XABAR  the  values  of  Z  ,  and  X'  calculated  by  the  heuristic 


//DAT* . INPBT  DC  • 
2 

200,225.30.10 
.05,.  10 
5,30 


STOCTKTIM  > C  F- FFOB  AE II IT y  FATTO  SAEPIJ9G  PIAF 


IAOL*  200.00  X1TPD=  225.00  SIGE1*  30.00 
AI PH  A  =  0. 05  BETA«0.10  KPBIBE*  10 

ZA  =  3.01  X1BAF=  214.26 


44 

44 

444 

444 

N 

4  * 

>44 

4 

444444444 

TL 

>44 

>44 

4 

TO 

>4444444444444444*  *4444 

♦•LOT  XE*F*PFOE  ACC* 

44 

444  1 

1 

44 

4 

123.57 

>44 

4 

301. 51 

>441 

44 

>4444441 

19  3.75 

>44 

4 

>4444444 

0.9e69 

444 

4 

2 

4 

348.26 

4 

506.77 

44 

195.21 

4 

0.9819 

4 

3 

4 

572.55 

4 

712.99 

44 

196.67 

4 

0.9752 

4 

u 

4 

796.63 

4 

917.21 

44 

198.12 

4 

0.9662 

4 

5 

4 

1021 . 12 

4 

1 121. 94 

44 

199.58 

4 

0.9541 

4 

6 

4 

1295. 41 

4 

1325.66 

44 

201.04 

4 

0.9381 

4 

7 

4 

1U69.70 

4 

1529. 89 

*4 

202. 50 

4 

C. 9  173 

4 

8 

4 

1653.98 

4 

1734. 11 

44 

203.96 

4 

0.8907 

4 

Q 

4 

1918.27 

4 

1938. 34 

44 

205.  92 

4 

0.8572 

4 

•>o 

4 

2112.56 

4 

2H2.  56 

44 

206. 87 

4 

0.8161 

4 

•  * 

444 

>4* 

>44 

>44 

208.33 

4 

0.7669 

4 

44 

209.79 

4 

0.7099 

4 

44 

211.25 

4 

0. 6459 

4 

44 

212.71 

4 

0. 5766 

• 

44 

214.17 

4 

0.5045 

4 

44 

215. 62 

4 

0.4321 

4 

*4 

217. 0« 

4 

0.  3623 

4 

4  4 

2’8. 54 

* 

0. 2976 

4 

4  4 

220.00 

4 

0. 2396 

4 

4  4 

221. 46 

4 

C.  1695 

4 

4* 

222. 92 

* 

0.  1474 

4 

4  4 

224.37 

4 

0.  1130 

4 

4  4 

225. 83 

4 

0.0656 

4 

44 

227.29 

4 

0.  064  1 

4 

4  4 

228.75 

4 

0.0476 

4 

4  4 

230. 21 

4 

0.0351 

4 

44 

231.67 

4 

0.0258 

4 

4  * 

233.  12 

4 

0. 0186 

• 

4  4 

234.5? 

4 

0.0137 

4 

4  * 

4  44 

4  4 

4  4 

44 

r«M 

236. 04 

• 

• 

C. 00°5 

• 

♦  •  4 

Figure  16.  Program  Output  Listing:  Format  2 


algorithm  and  used  as  arguments  in  the  calculation  of  the  SNPR  test 
plan.  The  column  headers  are  the  same  as  those  listed  for  the  SPR 
test  plan. 

The  following  are  the  data  cards  necessary  for  format  2  input: 

CARD  (2):  X ' AQL ,  X'LTPD,  a',  n' 

CARD  (3):  a,  8 

CARD  (4):  MAX,  NTERM 

Format  3:  SPR  and  SNPR  Test  Plans 

The  input  format  for  the  SPR  and  SNPR  test  plans  requires  speci¬ 
fication  of: 

(1)  x'AQL 

(2)  X'lTPD 

(3)  a* 

(4)  n' 

(5)  a 

(6)  8 

(7)  MAX 

(8)  NTRUNK 

(9)  NTERM 

Figure  17  shows  an  example  of  an  input  list  and  the  associated 
output.  In  the  output  listing,  a  header  which  echoes  common  and 
specific  input  for  each  test  plan.  X'  and  Z ,  are  listed  for  the  SNPR 
test  plan.  Column  headings  are  established  for  n  (sample  number), 
upper  and  lower  tolerance  limits  for  both  plans,  values  of  X'  for  the 
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//DAT  A . INPOT  Dt  » 
3 

233.22S.33, 13 
.05, . 1C 
5,20, 30 


*  • 


*  STCBEFTIAl  PRCEAEI1IT7  AND  NOJi-PFCBAEJLITI  PATIC  SASPING  PLANS  * 


A  A  Cl  =  200. OC  UTFD=  225.  00  SIGBA*  3C.0C 
SNFF  FLAN 

ALPRA=0 .05  PIT  A=  0. 13  NPEI«T=  10  2A*  3.01  X AB  AE=  214.26 

SPN  FIAN 

Al FHA  =0. 35  EITA=0. 10  TRUNCATION  OCCOBS  AT  N  =  20 


** 

00 

SNPRT 

000 

000 

00 

0000000000000000000 

SPB1 

*00 

00 

00*000+1 

i0 

*0 

00000001 

SNPPT 

t* 

* 

SPFT 

♦ 

* 

• 

V  »< 

TL 

0 

TO 

00 

TL 

* 

TO 

**LCT  XPAP*PPOE  ACC*PFOE  ACC* 

*» 

1  • 

122.97 

#** 

• 

00000000 

3DU.5U 

000 

00 

131.45 

*000 

0 

***#*#♦< 

316.55 

*00 

00 

00000001 

193.75 

0 

0 

*0000001 

0.987 

10 

0 

0.9e7 

2  • 

3ue.26 

• 

508.77 

00 

343.95 

0 

529.05 

00 

195.21 

0 

0.982 

0 

0.982 

3  • 

572.55 

• 

7  12.99 

00 

556.45 

0 

741.55 

00 

196.67 

0 

0.975 

0 

0.976 

u  • 

7S6.e3 

• 

917.21 

00 

768.95 

0 

954.05 

00 

198.  12 

0 

0.966 

0 

0.967 

5  • 

1021.  12 

• 

1121  .44 

00 

981 .45 

0 

1166.55 

00 

199.58 

* 

0.954 

0 

0.954 

6  * 

1245.41 

• 

1325.66 

00 

1193.95 

0 

1379.05 

00 

201.04 

0 

0.938 

0 

C.93P 

7  • 

1 46 $  . 70 

• 

1529.89 

00 

1 406.45 

0 

1591.55 

00 

202.50 

0 

0.917 

0 

0.  016 

8  * 

1693.98 

t 

1734. 11 

00 

1 618.95 

0 

1804. 05 

0  0 

203.96 

0 

0.891 

0 

o.epe 

9  * 

1 9 1  8 . 21 

• 

1938.34 

00 

1831.45 

0 

2016.55 

00 

205.42 

0 

0.857 

0 

0.e52 

10  • 

2142.56 

• 

2142.56 

00 

2043.95 

0 

2229.05 

00 

206.87 

0 

0.8  16 

0 

o.  eoe 

1 1  » 

2256.45 

0 

2441.55 

00 

208.33 

0 

0.767 

0 

0.764 

12  * 

00 

2468.95 

0 

2654. C5 

00 

209.79 

0 

0.710 

0 

0.693 

1?  * 

00 

2681  .45 

* 

2866. 55 

00 

211. 25 

0 

0.646 

0 

0.  624 

14  » 

00 

2893.95 

0 

3079. 05 

00 

212.71 

0 

0.577 

0 

0.552 

15  * 

00 

3106.45 

0 

3291. 55 

00 

214. 17 

0 

0.504 

0 

0.4-’7 

16  ♦ 

00 

3318.95 

0 

3504. 05 

00 

215.62 

0 

0.432 

0 

C.  405 

17  * 

00 

3531.45 

0 

3716. 55 

00 

2 1 08 

0 

0.362 

0 

0.  337 

18  * 

00 

3743.95 

0 

3920.05 

00 

218. 54 

0 

0.29e 

0 

0.27-> 

10  « 

00 

3956.45 

0 

u 14 1. 55 

*  0 

220.00 

0 

0.  240 

0 

0.224 

23  * 

00 

4168.95 

0 

4354. 05 

00 

221.46 

0 

0. '  89 

0 

o.  r? 

*  * 

*** 

00000000 

0  00 

* 000 

0  0 

222.92 

0 

0.  147 

0 

0.141 

00 

224. 37 

0 

0.  1  13 

0 

o.  m 

00 

225.83 

0 

0.386 

0 

3.0»7 

0  0 

227.29 

0 

0.364 

0 

O.Of  8 

00 

22e. 75 

0 

0.049 

0 

0.052 

0  0 

23C. 21 

0 

0.035 

0 

0.041 

00 

231.67 

0 

0.026 

0 

0.031 

0  * 

233.  12 

0 

3.319 

0 

3.  324 

0  0 

234. 58 

0 

0.0  14 

0 

0.019 

0  0 

00  * 

*0  0* 

0  0 

0  0 

236. C4 

0 

3.310 

0 

0.014 

Figure  17.  Program  Output  Listing: 
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OC  curves  and  the  associated  probabilities  of  acceptance  for  each  plan. 
The  following  are  the  data  cards  necessary  for  format  3  input: 

CARD  (2):  X'AQL>  X'LTPD.  o ' ,  n' 

CARD  (3):  a,  6 

CARD  (4):  MAX,  NTRUNK,  NTERM 

Program  Listing 

Following  is  a  listing  of  the  computer  code  which  generates  the 
test  plans.  It  consists  of  six  subroutines  and  a  driver  routine.  Each 
routine  is  commented  with  a  description  and  other  documentation. 
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***  20-0  (03/CU/R1--1BNU) 

C  H  A  IN  PRCGFA.F 

C  *•»•*»**««*♦ 

C  THIS  IS  THE  DRIVER  ROUTINE  FCF  TEE  FFCGBAil.  IT  El  ADS  DATA  AND  CALLS 
C  INDIVIDUAL  SU  EROUT  INIS  EASED  ON  A  PRESCRIBED  FCRHAT,  IPPINT. 

C  mHEFE  APE  THREE  POSSIELE  F0RBA7S  AS  DESCRIBED  IK  THF  OSEPS  GOICE. 

C  THE  SPF  TEST  PLAN  FOFEAT  CALLS  SUEFCOT3NE  SPFT  BHERI  TB£  TOLEEANCE 

C  LIHITS  FOP  THE  PLAN  ARE  CALCDLATED.  IT  TPTN  CALLS  SOBPCOTINE 

C  OUTPUT  BHEP.E  THE  CALCDIATICKS  AFI  COTFUT  ACCCKDING  TO  THE  SPECIFIED 

C  FOP  HAT .  THE  SNFF  TEST  FLAN  FCRHAT  CALLS  SUBROUTINE  SPRT  TO 

C  OBTAIN  THE  TOLERANCE  LIHITS  FCR  N=1  IN  THE  SPF  TEST  PLAN. 

C  NEXT,  SUBROUTINE  SSP  (SINGLE  SAHFLING  PLAN)  CALCULATES  THE 
C  SEED  VALDES  OF  ZA  ANC  XA  USED  El  THE  NEXT  CALLED  SUEROOTINE, 

C  SUPEOUTINE  SEARCH.  THIS  SUBROUTINE  RETURNS  THE  UPPER  AND  LCHFF 
C  TOLEEANCE  LIHITS  FCR  THE  SNFF  TEST  PLAN.  SUBROUTINE  OUTPUT  IS  THIN 
C  CALLED  TO  PRODUCE  THE  EECUIFEC  CUTFUT  LISTING.  THE  COHBINED  SPF 
C  AND  SNPF  TEST  PLAN  FORHAT  IS  THE  SAKE  AS  THAT  OF  PCFMAT  2  EXCEPT 

C  THE  OUTPUT  FORHAT  IS  ECDIFTED  TC  HANDLE  BOTH  TEST  FLAN  FORHATS. 

C 

C*********************************************************** ************ 

CCHHON  ALPHA, ASDH , B IT A , BSUH, S IGH A , XA , XAQI , XAS N PR , X ITPD , Z AS N FF , 

*  ZASFRT, IPFINT,K1X,NHAX,NPR1HE, NTEPH , NTR UN N , A S KCKV (100) 

*  ,OC (100) , CCS  N F F  (100) ,SNPPT1  (500) ,SNPRTU  (500) ,SPFTL(500)  , 

*  SPPTU (500) .XEFLCT  (100) 

FIAD,IPFINT 

GCTO (  1 , 2, V) , IPBINT 
C  PROCESSING  POFHAT  FCF  SEPT  ELAN. 

1  BEAD,XAOL,XLTPD,SIGHI 
BEAD, ALPHA, BETA 
BEAD, KTSUNK, NTEPH 
CALL  SPPT 

CALL  OUTPUT 
GOTO  5 

C  PPOCESSING  FOFHAT  POP  SNEFT  F l A  K . 

2  EFAE,XAOL,XITPD,SIGHI,NPEIHE 
BEAD, ALPHA, BETA 

PF AD, K AX, NTEPH 

C  ONLY  THE  TOLEEANCE  IIFITS  AT  N=1  ABE  NEEDED  FOB  THIS  FCBHAT. 

NTFUNK=  1 

3  CALL  SPPT 
CALL  SSP 
CALL  SEARCH 
CALL  SNPEOC 
CALL  OUTPUT 
GCTO  5 

C  PROCESSING  FORHAT  FCF  SPP1  AND  SNPFT  PLANS 
«  F.EAD.XAOl,  XLirC,SIGHA,NPFIHE 

FE  AC, ALFBA, BET  A 
PFAD, HAX, NTRUNF,  NTEFF 
GCTO  3 
5  STOP 

FNC 


SUBROUTINE  SSP 


C  *.***«•*.*«*** 

C  THIS  S  U3P OUT IN'E  IS  A  EEEEEF  ROUTINE  TO  SOBROOTI NE  SEARCH.  IT  CAICU- 
C  LATES  THE  SEED  VA10ES  CF  2AIFHI  OSJKG  TEF  VALDES  OF  TL(1)  AND  TUP) 

C  (SEE  CHAPTER  III  CF  TEESIS)  GENERATED  IK  SUBROUTINE  SPRT.  IT  ALSO 

C  CALCULATES  THE  SEE!  VALUE  OF  XBARPEIHE  TBC*  XABAE,  THE  INDIFFERENCE 

C  POINT  OP  ACCEPTANCE  IIKIT,  EOF  A  SINGLE  SAMPLING  PLAN. 

C 

C 4 

CCHSOK  ALPHA,  ASU2 , B IT A , E SUN , S IGH A , X A , XAQL , X ASN FE , I LTFD , Z ASNPF , 

*  ZASFRT,IFFINT,H AX  ,  N fl A  X ,  N PP IHE , KTERH . KTP ON K , A S NCR V (IDO) 

*  ,OC  (1 00) , CCS  NFS  (100) .SNFFTI  (500) .SNPRTU (500) ,SPETL (500)  , 

*  SPRTU (SCO) ,XERLC1  (100) 

C  CALCOLATE  THE  SEED  VAIUE  CF  2AIPHA,  ZASFFT. 

ZASPPT= (SPRTU (1) -SPFTL  (1) )/2/SIGEA 
C  CALCULATE  THE  REQUIRED  SAHPLE  NUMBER  OF  THE  SINGLE  SAPFLING  PLAN 
C  ENSSP  AND  THE  VAIOE  CF  THE  NEXT  HIGHEST  INTEGER,  NSSP. 

FNSSP=  ( (Z?BE  (ALPHA)  ♦ZPEB  (BETA) ) * SIGH A/AB S ( XA QL-XLT PD) ) **2 
NSSP=RNSSP 

C  AFTER  DETERMINING  VHICH  VALOE  IS  GREATER,  XAOL  AND  XLTFD,  CALCOLATE 
C  THE  INDIFFERENCE  FCINT  SITH  N=PNSSP  AND  N=NSSF.  TAPE  THE  AVERAGE 
C  OF  THESE  TV 0  VAIUES,  XA ,  AND  USE  IT  AS  THE  SEED,  XEAPPF IHE. 

IF  (XAQL.GE. XLTFD)  GOIO  10 

X AEAR  =  (ZPRB { ALFH A ) *SIG f A /P NSSP* *0 . 5)  ♦XAQL 
XAHI=  (ZPEB  (ALPEA) *S IGH A/ NS SP** 0 . 5) +XAQ1 
XA~AHAX1 (XABAP.XAHI) -AES (X  AB  AR-XABI) /2 
GOTO  20 

10  X A EAF= (-ZPRB(ALPBA)  *S I GH A/FN SSP** 0 . 5) ♦ XAQL 
XAHI*( -ZPRB (ALPHA) *SIGPA/NSSP**0. 5) ♦  XAQL 
X A*A H AX  1  (XABAP,XABI)-AES  (X AB AR- XA HI ) /2 
RETURN 
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subroutine  sprt 

c  **•*•••••••**•• 

C  THIS  F  CUT  INI  USES  THE  EQUATIONS  CFVELCFED  BY  PAID  (SEE  CHAPTER  II 
C  OF  THESIS)  FOP  THE  SPR  TEST.  IT  GENERATES  A  TEST  PLAN  FOP  THE 
C  SPECIFIED  INPUT  VlinES  CF  XAQI ,  XITFC.SIGBA,  ALPHA  ANt  FETA.  IT 

C  TRUNCATES  THE  PLAN  AT  TEE  INPUT  VALUE  OF  NTRT7  NF  CR  THE  DEFAULT 

C  VALUE  SBICH  IS  TV1CE  THE  VALUE  OF  THE  PAXIBOK  ASR  AT  IACL  CF 

C  XLTPD .  VALUES  FOE  THE  CPEPATING  C R A E ACT E FISTI C  (OC)  ANE  AVERAGE 

C  SABPLR  NUMBER  (ASN)  CORVES  ARE  ALSO  GFNFRATED. 

C 

CCMBON  ALPHA, A SUR , P ET A  ,  E SU B , S I G !> A , X  A , X  A CL, X  A £N PR, AL1P C , ZASNP P , 

*  ZASPRT, IPRINT  ,BAX  ,NBAX, NFEIEE, KTEFM, STRUNK,  ASNCFV (1  DO) 

*  ,OC(100) .OCSNPP  (100) , SNPBTL  (500) .SNPBTO  (500) , SFRTL (500) , 

*  SPRTU  (500)  .XERLCT  (100) 

C  DECLARATION  OF  FUNCTION  STATERINTS. 

ACC  (N) =H0*N*S; RSJ  (N) *H  1*A*S 
H  (THETA) = (XAQL«ILTPE-2*T FETA)/ (XLTFD-I ACL) 

PACC (TBETA)=  (A**H  (TPETA)  -1  )/  (A**B  (THETA)  -B**H  (TBETA)  ) 

ASN  (THETA)  =(H1*FACC  (THETA)  »(B0-B1)  )/(TEETA-S) 

C  CALCULATION  OF  EOUATICN  FAF1FETEFS. 

A=  (1-BETA) /ALPHA 
E=BETA/  (1-ALPHA) 

S= (XAQL+XLTPE) /2 
H0=SIGBA**2*ALCG  (E)/(XITED-IACL) 

H1=SIGBA**2*ALCG  (1) /  (X LTFD-I AOL) 

C  DETEKHINATION  OF  TFONCATICK  VALUE. 

IT  ERB  =  NTBUNK 

IF  (NTEUNK.LT. 0.5)  ITEBB=2*AHAX1  (ASN(XAQL)  ,ASN(XLTFD)) 
NTPUNK=ITEBB 

C  CALCULATION  OF  THE  TCIEFINCE  UB1TS. 

DC  130  1=1, ITE  RH 
SPFTL  (I )  =  ACC  (I) 

100  SPRTU(I)=PEJ(I) 

C  CALCULATION  OF  THE  STEP  FACTOR  USFD  F CF  THE  XBAHPFIHF  OF  TPE  OC 
C  ARC  ASN  CURVES.  TEE  1IK1TS  FOP  XBAPPBIHE  EXTEND  BEICND  XAOL  ANC 
C  X IT  PD  BY  AN  AMOUNT  DETERMINED  BY  THE  VALUES  OF  ALPHA  ANC  BETA. 

F ACT  1  =  ALPHA* 5; TACT2  =  EETA»5  ; FACTOR"  ( 1  +FACT 1 +FACT2) /KTEFF 
RINCR"=0;XLTPD1=XAQl-FICT1 »AES  (XACI-XLTPD)  ; NBAX=0 
C  THIS  SFGBENT  FILLS  TEE  ARRAY  XEAFLCT  WITF  THE  VALDES  OF  TEE  ARG0"ENTS 
C  FCP  CALCULATION  OF  TRI  CC  AND  ASN  CCRVES.  ARRAYS  CC  ANE  ASN  AFE 

C  FILLED  WITH  THE  ASSOCIATED  VALUES  CE  THE  PROBABILITY  CE  ACCEPTANCE 

C  ANC  ASN. 

DO  200  1=1 , NTEFH 

XBFLCT ( I) =XITPC1 ♦EINCFH 
OC  (I)=?ACC(XLTFD1*FINCRN) 

ASNCPV(I)  =  A  SK (YLTPC  1*FINCFB) 

C  CALCULATION  OF  THE  M  A  XI  F  OB  ASN. 

IF (ASSCRV (I) . GT. NBAX)  N BA X  =  A SNC P V  (I ) ♦ 1 
RINCPE=FINCF.*  +  FACTCF*  AES  (XACl-XITFC) 

RETURN 
END 


200 


cd  n  on  lp  cow  k>o-*  o  non  ooooo 


SUBROUTINE  SNPRT  (7 41 FH * , V BA R) 


THIS  S U3R  0UT1 NE  OSES  THE  IFFUT  VALUES  CE  XACL,  ZLTPE,  S1GHA  AND  NPRINE 
,  ALSO  THE  CALLING  ARGUMENTS  ZALPHA  »NE  XBAR  TC  GBRTRAT  E  AN  SNFP 
TEST  PLAN  WITH  ASSCCIATEE  ALPHA (ASOf)  ANE  EETA(ESUH). 
*««•»**«********••***»**•*»•****»*******»****»****»**»****•*»*»»**»••»» 
CCHKON  ALPHA, AS09 , E ETA , B SUR, SIGH A , XA , XAfil , X A SKFR , 1ITFD ,Z AS  HE F , 

*  ZASPRT, IFFINT  »K  AX , K  H AX  «  NPRIBE , KTEPH , RTF  0NN , AS  NCF V (133) 

*  ,OC  (100)  ,  CCS  A  F  F  (100)  .SNFFT1  (500)  .SN'PPTU  (500)  ,  SFETL  (500)  , 

*  SPRTO  (500)  .XERLCT  (100) 

CALCOLATE  THE  STEF  INTERVAL  (SEE  CHAPTER  III  CF  TRESIS). 

ENTRVL  =  ZALPHA*SIG»A/  (NFFIflE-1) 

CALCOLATE  THE  CONTROL  LIBITS  ANE  THE  ALPHA  AND  BETA  CF  THE  TEST  ANE 
STCPE  THE  RESOLTS  IN  THE  AFFAES  SNPRTD  AND  SNFPTL. 

DO  12  J=1 , NPFIEE 

RJ=J 

SNPFTO (J)  =XEAF*J+  (  (N FRIEE-J) »EKTFVL) 

SNPRTL  (J)  =  IEAF*J-  (  (NFFIHE-J) *ENTFVl) 

CALCULATE  ALPHA  USING  XAOL. 

DEVIAT= (SNPETO(J) - J* XA CL ) / (SOFT (RJ) * SIGH A) 

IF  (DEVIAT . GT . 0. 0 ) GCTC  1 
ALFH= 1-PRBZ (-EEVIAT) 

GOTC  2 

ALFH=PRBZ(EEVIAT) 

CALCOLATF  BETA  USING  XACI. 

DEV IBT= (SNFFTL(J) -J* X ACL) / (SOFT (RJ) * SIGH A) 

IF  (DEVIBT.GT. 0.0) GOTO  3 
PFT=PPBZ  (-EEVIBT) 

GOTO  4 

BET- 1 -PR  BZ  (DIVIBT) 

CALCULATE  ALPHA  SDR. 

CONTINUE 

IF  (J.GT.1.5JGOTC  5 
ASU  B= ALPH 
PRCD  =  1-  A  LPH-BET 
GOTO  6 

A5UH=ASUH*ALFH*EFCE 
PFCD=  PRO D*  (1-ALPE-EET) 

CALCULATE  ALPHA  USING  XLTPD. 

DEVIAT= (SNFRTU (J) - J*XLTP E) /  (S CFT  ( F J) *S IGH A) 

IP (DEVIAT. GT.O. 0) GCTC  7 
ALFH=1 -PRBZ  (-EEVIAT) 

GOTC  3 

A  LPH  =  PFB7  (DEVIAT) 

CALCULATE  EFT  A  USING  XLTPC. 

DEVIBT= (SNFFTI  ( J ) -0 *XLTP t ) /  (S g FT ( F J) *S IGH A) 

IF (DEVIBT. GT . C. C) GCTC  9 
BET=P?BZ  (-CEVIET) 

GOTO  10 

9  BET=1-PP PZ (DEVIBT) 

C  CALCULATE  BETA  SUE. 

10  IF  (J.  GT.  1.  5)  GCTC  1 1 
BSUM=  BET 

PFCDB=1-ALFH-EET 
GCTC  12 

11  E5U*=?SUR*FET*PRCEE 

PR  C3B=  FF  CP  E*  (1-A1FH-EET) 

CONTINUE 
RETURN 
*NE 
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SUBROUTINE  SNPFOC 

•**•***»**••**••• 

THIS  SUBROUTINE  CAICULATES  VALUES  fCF  THE  CC  CURVE  OP  AM  SNPB  TEST 
PEAK. 


CCNNON  ALPHA,  ASD  S  ,  B  ET  A  ,  ESI!  E ,  SIGH  A  ,  A  A  ,  X  ACL  ,  X  A  SN  PP  ,  II  TP  T  ,  ZA  S  K  FF  , 

*  ZASPRT,3FFlNT.HAX,NnAX.NPFIFE,KTEFN, NTRUNK.ASNCBV (IOC) 

*  , OC (100)  , 0  CS  NPF  (100) , SKPB TL  (500)  ,SNPRTU  (500) ,SFRTL(500)  , 

*  SPRTU  (500) ,XEFLCI (100) 

CALCtTLATE  VALUES  FOB  T  BE  CC  CURVE  0SIKG  THE  AREA?  XBPJOT(IOO). 

DC  6  K=1,NTERH 

CALCULATE  THE  FBOEAEILITT  OF  ACCEPTANCE  FOR  A  SFECIFIC  XPPLOT. 

DO  7  J=1,NPSIHF 
HJ=J 

CALCULATE  ALPHA  USING  XERLOT  (K)  . 

DEVIAT=  (SNFRT0 (J) -0*XERLCT (K) )/ (SCRT (RJ)  *SIGMA) 

IF  (CEVIAT. GT. 0.  0)GCTC  2 
ALFH=1-PRBZ(-DEVIAT) 

GOTO  3 

ALFfi=PBBZ(EFVIAT) 

CALCULATE  BETA  USING  XEELCT(K). 

DEVIBT= (SNPPTL(J) -J*XBBLCT (X) ) / (SQPT (RJ) *STGSA) 

IE (DEVIBT. GT. 0. 0) GC1C  K 
BET=PEBZ (-EEVIBT) 

GOTO  5 

BET=1-PRBZ  (EFVIET) 

CALCULATE  BET  SUH. 

IP  (J.  GT.  1.  5)  GOTO  6 
BETS08=BET 
PFODB= 1- ALFB-EET 
GOTO  7 

6  BFTSUN=EETSDH*BFT*PFCDE 
PFCDB=FBODE*  (1- ALFB-EET) 

7  CONTINUE 

8  CCSNPP  ( K ) =E ETS OH 
RETURN 

FKE 
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SUBROUTINE  SEARCH 

C  •**«•*•••**•*•»*• 

C  THIS  SUBROUTINE  IN  CONJUNCTION  WITH  TBE  CTHEP  SUBROUTINES  BINDS  THE 
C  PARAMETERS  ZAIPHA  AND  XBARPRIHE  WHICH  POR  A  SPECIFIC  VALUE  OF 
C  NPRIHE  YIELD  AN  SNPR  TEST  PLAN  WITH  ASOH=ALPRA  AND  BSUN=BETA. 

C  THIS  IS  DONE  USING  THE  HEURISTIC  ALGCFITRH  DEVELOPED  IN  CHAPTER 
C  III.  USING  THE  SEED  VALDES  OP  ZAIPHA  AND  XBARPRIHE  FROM  SOERCDTINE 
C  SSP,  A  MODIFIED  VALUE  OF  IBARPRIHE  (XBARPRIHE  (NPRIHE) )  IS  PCOND. 

C  AFTER  A  CHECK  TO  SEE  IF  THE  VALUES  CF  ASUH  AND  ESUH  ARE  WITHIN  THE 
C  DESIPED  TOLERANCE,  THE  SUEPOUTINE  IS  TERMINATED  (TCIEPANCE  MET)  CR 

C  THE  VALUES  CF  ZALPHA  AND  XEARFRIME  IRE  MODIFIED  (TOLERANCE  NOT 

C  MET).  A  RECALCULATION  OF  ASUS  ARC  BSD 3  WITH  THE  MODIFIED  CALLING 
C  ARGUMENTS  IS  FOLLOWED  EY  ANOTEER  TOLERANCE  CHECK.  THE  FROCESS 

C  CONTINUES  UNTIL  THE  VALUES  CP  ASUH  AND  ESUH  ARE  WITHIN  TOLERANCE  CR 

C  THE  MAXIMUM  NDMEEF  OF  ITERATIONS  (MAX)  HAVE  BEEN  PEBFOFHED. 

C 

C  **************************** ••*•*****♦**• ****************************** 

CC.M30N  ALPHA,  ASUH,  BETA,  BSUH,  SIGMA,  X A , X  A QL, X  ASH  PR .  XLTPD , ZASNPB , 

«  ZASPRT, IPR1NT,HAX,NHAX,NPRI3E,NTERH, NTBUN K , AS RC RV ( 1 00) 

*  ,OC  (100)  .CCS  HEP  ( 100)  ,SNEETL  (500)  ,SNPRTU (500) ,SPRTL (500) , 

*  SPRTU  (500)  .IEBLCT  (100) 

R N FRI M= NPRIHE ; RN MAX =NH AX 

C  CALCULATE  Z ALPHA  (N FRI ME ) 

Z  AT=Z  ASPRT* (1.0*  (1.0-BNPEIB/RRHAX) ) ; IBBTHP=X A ; BTIHE= 0 
C  CHECK  TO  SEE  IF  ASUH  AND  ESUH  ARE  WITHIN  0.001  CF  THE  SPECIFIED 
C  ALPHA  AND  BETA  AFTER  CALLING  SUBROUTINE  SNPRT. 

1  CALL  SSPRT (Z AT, X BBT HP) 

IF  (ABS (ALFHA-ASUH) . LE. 0 . 00  1. AND. ABS  (BETA-BSUH)  .LE.O. 00  1) GOTO  50 
IF  (NTIHE.GT. MAX)  GO  TO  50 

C  CALCULATE  THE  DISTANCES  OF  ASOH  AND  ESUH  PRCH  ALPHA  AND  BETA. 

Dlf-ABS ( (ASUH- ALPHA)  -  (ESOH-BETA) ) 

C  DETERMINE  IF  ASUH  HOST  INCREASE  OR  DECREASE  AND  ALTER  ACCORDINGLY. 

C  CALCULATE  NEW  ASUM  AND  ESUH. 

IF  ((ASOH-ALPHA) , LT.  (BSOH-B ETA) ) GOTO  10 
SIGN* 1 
GOTO  20 
10  SIGN*- 1 

0  AMARK  =  DIF* ( 1 / ( BET A/ A LPB A  ♦  1 ) ) 

ALTER  XBARPRIHE  EY  A  FIXED  AHOONT. 

SENEST=SIG MA/1 00.0 

XBRTH 1  =  XBRTHP* SEN  ESI* SIGN; ASUHPV  =  ASDH 
CALCULATE  THE  EFFECT  CF  THE  ALTEFATIC*. 

CALL  SNPRT  (Z AT, XERTH1) 

CALCULATE  THE  PROPORTIONATE  VALUE  BY  WBICH  XBARPBIHE  MUST  BE  ALTERED 
TO  MEET  THE  OBJECTIVE  VALUE  OF  ALFBJ. 

SENCAL  =  SENEST* AM  A RK / A B S  (  A S Ufl P V- A S UH) 

XEPTMP*XBPTMF»  SEKCAI ‘SIGN 
CALL  SNPRT  (Z AT,  XESTHP) 

CHECK  TO  SEE  IF  ASUM  AND  ESUH  ARE  SIT1N  DIBITS. 

I?  (ABS  (ALPHA -A  SUM).  LE  .  0 . 00  1 .  AN  D.  AES  (EETA-8SUH)  .  LE .  0 . 00  1)  GOTO  50 
DETERMINE  IF  ASUH  MOST  INCREASE  OB  DECREASF.  ALTER  ZAT  ACCORDINGLY 
FOR  A  NEW  ASOH  AND  BSUH. 

IF  (  (ASOM-ALPHA) . IT. 0.0)GCTC  30 
SIGN*  1 
GOTO  UO 
30  SIGN*-  1 

1*0  AM  ARK*  AES  (  ASUM-ALPHA) 

C  ALTER  ZALPHA  BY  A  FIXED  AMOUNT. 

ZATT=Z  AT  *0.0  1*  SIGN;  ASUNEV=ASCH 
CALL  SNPRT  (ZATT, XERTMP) 


U  U  U  u  U  U  U  *-  *-  u  <N 
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SOERO  OT I N£  OCTFOT 

THIS  SUBROUTINE  FORMATS  AN  I  COTPOTS  THE  RES0L1S  ACCORDING  TC  THE 
REQUESTED  FORMAT. 


CCKMON  ALPH A, ASU  “  ,  B  ITA  ,  B  SO  M, SIG H A , I A , I AOL , X A SN PR , I ITTD , ZAS NF F , 

*  ZASPBT, IPFINT.MAX.NHA X.NPPIBE, KTER M , NTR ONK , ASNCFV (10D) 

*  ,OC  (100) ,CCSNFF  (100) , SNFFTL  (500) .SNPPTU (500), SPPTL (500) , 

*  SPSTO  (SCO) ,XERICT (100) 

DETERMINE  THE  OOTFOT  FCFHAT. 

ITERM=NTRONK 
IF  (IPKINT.  EQ.  1 )  C-OTO  1 
IF  (IPRINT. EQ. 2) GOTO  2 
IF  (IPFIKT.  EQ.  3)GCTC  3 
FOPMAT  FOR  THE  SPF  Ft A K . 

WRITE ( 6, IDO) 

WRITE  (6,101) 

K R IT E  ( 6 ,  102) 

WRITE  (6,101) 

WRITE  (6,  100) 

PRINT,'  « 

FEINT,'  • 

WRITE  (6,  103) 

PRINT,'  • 

WRITE  (6, 10tt)  XACl.XLTPD, SIGMA 
W  R  IT  E  ( 6  ,  1 3  5)  ALPHA, EFTA, 1TER8 
PRINT,'  • 

WRITE  (6,  10  3) 

WRITE (6, 103) 

PRINT,'  • 

WR IT E  ( 6 ,  106) 

PRINT,'  • 

WRITE  (6,103) 

DO  10  1=  1 ,  IT  FR  S 

0  WRITE (6,1 07)  I  ,SFFTI  (I) ,SPFTU  (I)  ,XEEIC7 (I)  ,OC  (I) . ASNCFV (I) 

WRITE  (6, 108)  XFRICT  (ITFRR+1)  ,OC  ( ITFR  B ♦ 1 )  ,A SNCP V (ITERH* 1) 

ITEMP  1=ITERM*2 
DO  11  I=ITEMP1 .NTERF 

1  WPITE  (  6,  109)  X  EPLOT  (I)  , OC  ( 1)  ,  A S NCR X  ( I) 

KFITE (6,103) 

GOTO  1000 

FOPMAT  FOP  THE  SN  FF  FLAN. 

WRITE (6, 200) 

WRITE  (6,201) 

WFITE  (6,202) 

WPITE  (6,201) 

WPITE (6,230) 

PRINT.'  • 

PRINT,'  ' 

WRITE (6,200) 

PRINT,'  ' 

WRITE (6, ID  U)  XACL.FLTPD, SIGM A 
WRITE (6, 20  U)  ASOM, BSOM.NPRIRE 
WPITE (6,203)  2ASNFP.XASNFF 
FRIST,'  ' 

WPITE  (6,200) 

VPITF  (6,200) 

PRINT,'  ' 

WRITE  (6,  2D  5) 
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PRINT,'  • 

WRIT!  (6,200) 

DC  20  1= 1, NPRI f E 

20  WF ITS  (6,206)  I , SNFRT 1  (I )  ,S KPFTO  (I) ,X ERI OT (I)  , CCS NPF (I) 

WRITE (6, 207)  X  ERECT (KPBIKE«1) .OCSNPP  (NPRIBE*1) 

ITEBP1«NPRTflE*2 

DC  21  I=ITEBPi ,NTEBS 

21  WE ITE  (6,20  8)  XEFICT  (I)  .CCSNPF (I) 

NR  IT  E  (6, 200) 

GOTO  1CC0 

C  POFBAT  FOB  THE  SNFF  »ND  SIB  P1ABS. 

3  V RITE  (6, 300) 

WRITE  (6,301) 

WRITE  (6,302) 

WRITE  ( 6 , 30  1) 

WPITE (6,300) 

PRINT,'  • 

PRINT,'  • 

WRITE  (6, 30  3) 

PRINT,'  • 

WPITE  (6,304)  TICI,X1TFC,SIGB* 

WRITE  (6, 30  5) 

WRITE  (6,306)  ASOF,BSOE,NFRIBI,ZASNFR,IASNPR 
WRITE  (6,307) 

WPITE  (6,308)  A1PBA, BETA, ITERS 
PRINT,'  • 

WRITE  (6,30  3) 

WRITE (6,303) 

PRINT,'  • 

WRITE  (6,  3091) 

WPIT  E (6, 309) 

PRINT,'  • 

WPITE  (6,303) 

DO  30  I*  1, NPRIBE 

30  WRITE  (6,310)  I  ,SRPRTI  (I)  ,SNPFTO  (I)  ,SFPTL (I) , SPPTD  (I) ,T  ERLOT  (I)  ,OCS 
•PPF (I) ,CC(I) 

I=NPRI HE* 1 

WF ITE (6,311)  I ,  SFPTI  (I) ,SPFT0  (I)  ,XEPICT(I)  ,OCSNPR(I) ,OC(I) 

ITEBP=RFRIHE*2 

DO  31  I=ITEBP,ITEFfi 

31  WRITE (6, 312)  I,SFPTI(I) ,SPRTD(I)  , IERIOT  (I)  ,CCSNPP(I)  ,CC(I) 

WRITE  (6,313)  XERICT  (ITEEF*  1)  ,CCSKFF  (ITEPE*1) ,OC(ITEPB*1) 
ITE(1P1=ITERM*2 

DO  32  I=IT  EBP  1 , NT  ER  B 

32  WRITE  (6,314)  XEFICT  (I)  ,CCSSPF  (I)  ,CC  (I) 

WPITE  (6,  30 3) 

GOTO  1000 

100  FCFflAT  (6X,  46  ('  ••)  ) 

101  FORBAT (6X, '*' ,  44  X ,' •*) 

102  FCRBAT  (6X,  '  *  SEQUENTIAL  FFOBABIIITY  RATIO  SABPUNG  PLAN  *•) 

103  FOBBAT<1X, 56  (•••)) 

104  FORBAT  (6X,'XA0I=’,P7.2,2X,  'X  IT  FB=  ’  ,  E7 . 2, 2X,  •  SIGH  .»  =  •  ,  16 . 2) 

105  FCRBAT(6X,'ALPPA=',E4.2,2X,*PETA«',F4.2,3X,'TPnNCATE  AT  X=  *,I4) 

106  FORBAT  (IX, ••', 2X,'N' ,2X,  ••', WX, Ml', flX,'  T0',4X,'*ICT  XBAB*FRCE  AC 
*C*  ASN  *•) 

107  FORBAT  (1 X, • *  * , 14 , IX , • * • . F9 .2 , •  *',F9.2,'  *',17.2,'  *  ',F6.4,'  *• , 

*E6. 2,  •  ••) 

108  FOPB AT  ( IX, 29  ( •  •• )  ,F7.2,'  «  »,F€,4,»  »',F6.2,'  *•) 

ICO  FCRBAT  (IX, •*• ,26X,' •••  ,F7.2,'  *  ',16.4,'  *»,F6.2,'  *•) 

200  FCP.B  AT  (  IX  ,  50  ( '  •  ' )  ) 


FCRHAT (  IX, •*», USX  , ••  •) 

FORSAT  (IX, ' *  SECUENTIAl  ECS-F FCBAEIIITI  FATIC  SASPLING  PLAN  *») 

FC  SB  AT  (6T, 'ZA=',F6.2,2X,  'XAEAF=',F7. 2) 

FCFSAT (6X, • ALPEA*', F0.2.2X,'  EETA=*  ,F4. 2,3X, • KPRI SE= ' , I U) 

FOPB  AT  (IX,  *  *  N  *',«X,'TL',«X,'*' ,N7,'T0' ,«X,'**ICT  XBAR*  FFCP  ACC 
*•') 

FCFSAT  (IX, •*', J9,'  *  '  ,  F9  .2  ,  *  *',T9.2,'  **',F7.2,'  •  *,F6.U,'  *•> 
F0EBA7  (IX. 30  (•  *•)  ,F7.2,  •  *  ',F6.ll,'  »•) 
FOFSAT(1X,'*',27l,'««',F7.2,  •  •  ',F6.9,'  *•) 

FCESAT  (8X, 68  (•♦')) 

FOPS  AT  (8X,  *••,€«,  •  •  •) 

FORSAT  (8X,‘*  SECDEFTIA1  PPOEAEIIITT  AND  NON-PRO EAEILTTT  RATIO  SAS 
•FING  FLANS  *•) 

FORSAT (1X,82(*  «') ) 

FORSAT  (II, 'XAQL=',F7.2,2X, *X tTPC= • , f 7. 2,21 , • SIGH A= • , F6. 2) 

FCFRAT  ( 10X  , *  SNPH  PLAN') 

FOBS  AT  (IX, • ALPBA=', !«. 2, 2X, « EET A* • , F« . 2, 2X , • HPEIHE* • ,1 U ,3X ,  •  ZA= • , E 
*6.2,2X,'XABAP=',F7.2) 

FCFBAT  (10X,'SPF  PLAN') 

FORSAT (IX, 'ALPHA*', FA. 2, 2X, *EETA=  '  , E« . 2, 3X, • TRUNCATION  OCCOPS  AT  N 
•»',I«) 

FCRSAT (17X, 'SNPRT*, 8X, ••*',9X, 'SPET*  ,20X , • SNPRT' ,2X, • *  SPRT') 
FORSAT  (IX.  •*  A  »»•  ,«X,  'TI*  ,41,  'TO  ',  UI,  •**',  «X,'TL',  «X.  •  •• 

•,UX,'TO',HX,'**LCT  IE  A  E*  FRCB  ACC*PFCE  ACC*') 

FCFSAT  (IX, '*', I#, '  ,F9. 2,  '  *',F9.2,'  **',F9.2,‘  *',F9.2,'  **',F 

*7.2,'  *  • ,  F5.  3  ,  •  •  ♦,  F5 .3,  •  *•) 

FCBSAT  (IX, *•' ,T4,1X  ,25  ( • *  •) , F9 . 2,  •  *',F9.2,»  ♦*',F7.2,'  *  *,F5.3, 
*  *,F5.  3,  • 

?ORSAT(7X,'*',I«,'  •••  ,21K,'**',I9.2.'  *',F9.2,»  **',T7.2,'  *  ',15 
*.3,'  *  • , F5 . 3  , '  *•) 

FCBSAT  (IX,  5H  ( •  *')  ,  F7 . 2,  '  •  ',FS.3,'  •  '^5.3,'  *•) 

FORSAT  (1X,'*',51X,'**<  , F7 . 2 ,  '  *  ».I5.3,'  *  *,F5.3,'  *•) 

RETURN 


APPENDIX  B 


SELECT  SNPR  TEST  PLANS 

The  following  pages  contain  SNPR  Test  Plans  for  twelve  sets  of 
sample  data.  In  each  case  X 1  Aqj_  =  45  and  =  60.  Three  different 

values  are  used  for  a';  15,  22.5  and  30  (these  values  represent  1,  1.5 
and  2  times  the  difference  ( ^ ' AQL  "  *'lTPD^‘  fr°r  eac*1  va^ue  two 

plans  are  generated  with  n'  equal  to  the  maximum  ASN  of  the  corresponding 
SPR  test  and  approximately  70%  of  that  value.  Two  sets  of  plans  for 
a  =  .05,  6  =  .10  and  a  =  .01,  8  =  .01  are  generated  in  each  of  the 


above  cases. 
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e  v  ur  c  o  c  o  c  c  r*  a 

c  a  «_<  l»  ‘  «-'  r»  __  w  <MCjr-0  **  P-  f><c:u  ~,0w*0,""p',,'±-«  r '  x  P 

iT  ^  rJ  ff  w  r  *■  ;  /  r»  r  p*  3  r\  c  £  j  ^  p  ai  f”  •”  cr  r"'*-'  3  i*  tN  »■  v  a.  r*-  v*.  u’u”^  r  f*i  rri  r.  ^  <\  r-  *-  *— 

r-r-r-^^siw^aw  u  ;y  cj  ^  f  /'r‘«*»mtNr\.rNiNfVT-*-*—  ^-»-^«-c'0<-;ec:oc©c:r3eOCC5o*xcO'~'T 


C  C  f  ^  C  »j  c  o  C  CC  I.C  c  C  r)  c  u  c  •  J  O  c  c  C  o  C  ‘  ■>  ^  c  o  u  *  ’  C  ^  C  C  O  c  ^  C-' 
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